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As primitives for entanglement generation, controlled phase gates take a central role in quantum computing. 
Especially in ideas realizing instances of quantum computation in linear optical gate arrays a closer look can be 
rewarding. In such architectures, all effective non-linearities are induced by measurements: Hence the probabil- 
ity of success is a crucial parameter of such quantum gates. In this note, we discuss this question for controlled 
phase gates that implement an arbitrary phase with one and two control qubits. Within the class of post-selected 
gates in dual-rail encoding with vacuum ancillas we identify the optimal success probabilities. We construct 
networks that allow for an implementation by means of todays experimental capabilities in detail. The methods 
employed here appear specifically useful with the advent of integrated linear optical circuits, providing stable 
interferometers on monolithic structures. 



PACS numbers: 03.67.Lx,42.50.Ex,42.50.Dv 

I. INTRODUCTION 

Linear optical architectures offer the potential for reliable 
realizations of small-scale quantum computing IH In the re- 
cent past, numerous proof-of-principle demonstrations have 
relied on the precise state manipulation that is available us- 
ing linear optical elements. The development of better and 
brighter sources with good mode quality as well as new 
types of detectors have opened up new perspectives |2] in 
state preparation and manipulation, for six or more photons. 
Specifically, integrated optical circuits allow for state manip- 
ulation with little mode matching problems in interferome- 
ters H. 

Naturally, significant efforts has been devoted towards re- 
alizing instances of quantum gates. As primitives for such 
small-scale computing, two-partite quantum gates delivering 
a controlled phase-shift of (p have already been experimen- 
tally demonstrated (see Refs. \4i-^ for ip = tt and Ref. |01 for 
general phases). In this note, we focus on linear optical imple- 
mentations of phase gates with arbitrary phases. In particular, 
we will ask when arbitrary phases can be realised in the first 
place, and — one of the main figures of merit in linear optical 
applications — what the optimum probabilities of success are, 
as any non-linear map is necessarily probabilistic. A realisa- 
tion of such gates seems interesting from the perspective of 

(i) gaining an understanding of the probabilistic character 
of quantum gates as well as 

(ii) serving as a proof of principle realisation of a kind of 
quantum gate that has several applications in linear op- 
tical quantum information processing. 

As far as the first aspect is concerned, one may well ex- 
pect that there is a trade-off between the notorious problem 
of having a small probability of success and the phase that 
is being realised in the gate. In fact, the study in Ref. [8] 
suggests exactly such a behavior: the presented upper bounds 
to the probability of success increase from the minimum at 



Ps{p = tt) = 1/4 to Ps{0) = 1. To investigate such a trade- 
off is interesting in its own right and helps in building intuition 
concerning the probabilistic behavior of linear optics. One 
may well develop the intuition that "large phases are costly" 
as far as the probability of success and hence the overhead or 
repetition are concerned. 

Further, concerning the second aspect, there are several ap- 
plications for which such a trade-off is relevant. In linear op- 
tical architectures, it may be a good idea to have a smaller 
phase, if one only has higher success probabilities. The new 
measurement-based quantum computational models for 
example offer this perspective: One does not have to have con- 
trolled TT phase gates to prepare cluster states, but one would 
in principle also get away with smaller phases. This may well 
(but does not have to be) a significant advantage when prepar- 
ing resources for measurement-based quantum computing dif- 
ferent from cluster states [l^ . 

Of course, in standard gate-based quantum computing, one 
will typically encounter all kinds of controlled phase gates. 
For example, in the quantum Fourier transform ifllfl . one has 
to implement several controlled phase gates. They can again 
be decomposed into other sets of universal gates (like CNOT 
or CZ and local unitaries). But, in terms of resource require- 
ments, it is obviously an advantage to directly implement the 
relevant quantum gates with phases in the range < p < it. 
There are also interesting trade-offs between resource require- 
ments and success probabilities in a number of related con- 
texts, like non-local gates in distributed quantum computa- 
tion lfl2ifl3ll . Refs. ifisll . for example, study distributed con- 
trolled phase-gates which would need less entanglement and 
succeed with a higher probability. In the field of linear op- 
tics gates, numerical results on direct implementation of arbi- 
trary two-qubit gates are known (see Ref. ||T4}| and references 
therein). 

Instead of resorting to decompositions in the circuit model, 
one could gain from implementing unitaries in a fashion "nat- 
ural" to the respective architecture at hand. In the case of 
linear optics it means to leave the computational sub-space 
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given by the encoding of the qubits for the sake of taking a 
"shortcut" through higher dimensions 101 ■ Given that the fun- 
damental information carriers are implemented using bosonic 
modes, this will occur when mixing those modes in beam 
splitter networks and is an inherent feature of genuine lin- 
ear optics implementations, in contrast to decompositions into 
standard gate sets. 

Here we will study post-selected gates, so not genuine 
"event-ready" quantum gates, but — as is common in linear op- 
tical architectures at least to date — those where one measures 
the output modes and whether the gate actually succeeded 
is determined only a posteriori by accepting only those out- 
comes which lie in the computational dual-rail subspace of 
the Hilbert space of n photons on 2n modes. Incorporating 
less constraints, these gates concern only a smaller number of 
modes and are still within reach of current experiments. In 
principle non-demolition measurements of the output would 
be required for an event-ready gate. 



respectively. Since we are restricted to n < 2 modes, these 
four quantities determine the action of the core completely. 
With the constraint 

1 - 2siii2(t?) = sin2(i5) (7) 

that ensures equal single- and two-photon amplitudes (equal 
probabilities for all dual-rail states), only the two solutions 

■d = ± arcsin(3~^/^) are possible, giving rise to (p ^ 0,tt, 
respectively. 

Hence, one finds that in this way, one can implement quan- 
tum phase gates, but only two different ones: One is not doing 
anything, and the other ones effect is a controlled phase of tt. 
This is exactly the gate of Refs. flBt]. In other words, with- 
out invoking at least a single additional mode, one can not go 
beyond the known 7r-phase in this fashion. 



B. Arbitrary phases 



II. CONTROLLED PHASE GATES 
A. Single beam splitter 

In a post-selected phase gate on four modes in the standard 
dual-rail encoding, two of the modes are merely involved as 
"by-standers", in that their amplitude is compensated in ex- 
actly the same fashion as in Refs. 10-01 ■ In this section, we 
will hence concentrate on two modes forming the "core" of 
the scheme, giving rise to a two-qubit dual-rail phase gate on 
four physical modes. The core itself may be regarded as a 
single-rail phase gate in its own right. Later we will see that 
not breaking the network into a core and by-stander modes 
will not give any advantage. 

Similar to Ref. we will briefly investigate the conse- 
quences of simply having a single beam splitter forming the 
core of the quantum gate. The action on the photonic creation 
operators of the two involved modes it is mixing is described 
by the matrix 
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with 



However, we can extend this scheme: The restriction to 
unitary two-mode beam splitters can be relaxed. Instead of 
starting with U £ SU{2), we use an arbitrary matrix A G 
(1-12x2 'pjjgjj yyjjj embed the two-mode matrix into a higher 
dimensional unitary, such that an appropriately rescaled A 
forms a principle submatrix of a larger unitary matrix A'. The 
optimal rescaling is simply dictated by the largest singular 
value of A. We will see that in the two-mode case only a 
single additional mode is already the most general extension, 
so the full set-up would consist of a transformation on three 
modes involving at most three beam splitters. In this class 
of gates, for each (p, the one with the optimal probability of 
success Ps{(p) can be identified. 

Observation 1 (Optimal post-selected dual-rail controlled 
phase gate). Consider linear optics, an arbitrary number of 
auxiliary vacuum modes and photon number resolving detec- 
tors. When post-selecting the state of the signal modes onto 
the computational sub-space and the auxiliary modes onto the 
vacuum, the optimal network on four modes implementing the 
gate represented in the computational basis of two dual-rail 
qubits by U = diag(l, 1, 1, e*"^), (p G [0,7r], has a success 
probability ( shown in Fig. [7} of 
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sm{'d) cos{'d) 
- cos{'d) sin(i?) 
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and appropriate phases 0i,^2 G [0, 27r) and mixing angle 
i9 e [0, 2tt). The phases can also be realised deterministi- 
cally by local operations on the dual-rail qubits, which leaves 
the relevant part of the gate U' = B. The matrix elements 
of the unitary U' belonging to U' for vacuum, single photon 
operation, and the two-photon component read 



(0,0|Z^'|0,0) = 1, 
(1,0|Z^'|1,0) - =sin(z?), 
(0,l|Zi'|0,l) = ^2,2 - sin(i?), 
(1, l\U'\l, 1) = pcr(A) = cos(2i9) 



(3) 
(4) 
(5) 

2sm'^^, (6) 



Ps{p)= 1 + 2 
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Proof: In order to find ps — the same for all possible in- 
put states — we will first construct the linear transformation of 
the relevant creation operators and identify the optimal unitary 
extension afterwards. 

• Two-mode transformation: The two-mode transforma- 
tion resulting from solving the equations (ISll-© im- 
posed by the gate we want to build is 



A = v'J' 



y/x l/x 



(9) 
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FIG. 1. Optimal success probability Ps{'p) of phase gates with vac- 
uum ancillas (one vacuum ancilla is already optimal) vs. the phase 
ip (solid line). A\.ip = -k the resuk of Refs. HlJl], Ps(vr) = 1/9, is 
reproduced. The intuitive assumption of a monotonous Psif) is not 
fulfilled: indeed, the success probability is worse than 1/9 in the in- 
terval 7r/3 < 1/3 < TT. Due to phases ip < it not being implementable 
with a single beam splitter, the additional unitary extension requires 
further measurements and therefore decreases the probability of suc- 
cess near tp = -k. 



X and y are free non-zero complex parameters. By writ- 
ing 



^ = py^diag(a,a-i) 



1 e*"^ - 1 
1 1 



diag(6,6-i), (10) 



with a = a;?/~^/^ and b ~ y^l"^ we see that the singular 
values of A only depend on |ap and so not on the 
phases of a,6, and x and y. 

The general solution to the dual-rail problem is actu- 
ally composed of the transformation A together with 
appropriate damping of the by-standers: the probabil- 
ity of success can not be enhanced by considering a full 
transformation on all four modes. This can be seen by 
writing the polynomial system given by the dual-rail 
problem similar to (O-®, consisting of 16 quadratic 
equations in the matrix elements of B e (D*^"*. It turns 
out that by permuting modes and appropriate variable 
substitutions all solutions can be brought into the form 
B o^\2® A, consisting of the two-mode core given in 
Eqn. (|9]l and two by-passed modes. 

• Optimal extension: Given the 2x2 matrix A that re- 
alises the transformation we are looking for, the optimal 
unitary extensions can be identified. Let us extend the 
first and second row vectors (denoted by A^ and A-2) to 
dimension 3 by appending ^13 and ^2,3, respectively, 
in such a way as to allow for unitarity of the extended 
matrix, N e SU{i). 

To see why a dimension of three is already sufficient 
consider a linear transformation of the creation opera- 
tors of n modes, described by its a (not necessarily uni- 



tary) matrix A. By using the singular value decomposi- 
tion (SVD) it can be decomposed as A = F • • W^^ 
where V and W are unitary (and therefore have imme- 
diate interpretations as physical beam splitter matrices 
themselves), and D = {di, . . . , (i„} is a diagonal ma- 
trix with real non-negative entries di > d2 > ■ ■ ■ > dn, 
the singular values. In terms of linear optics D can be 
interpreted as mixing each mode k — 1, . . . , n with an 
additional mode n + fc in the vacuum state which will be 
post-selected in the vacuum afterwards lfl5i[l6ll . Then, 
dk describes the transmittivity of the beam splitter used 
to couple modes k and Uk- Without loss of general- 
ity one can assume di = 1, vacuum mixing for the 
first mode. This can be achieved by rescaling with the 
inverse of the largest singular value, so A 1— > d^^A, 
which implies d^ 1— > d^^dk- Note that such a "global" 
rescaling of A does not change the post-selected action 
on the computational sub-space but only the success 
probability of it according to ps i—> d^'^^ps- Therefore, 
in general, there are n — 1 additional vacuum modes re- 
quired to extend an n-mode linear transformation to a 
unitary, and thus physical, network. Also please note, 
that the constraint di = 1 has to be taken into account 
for any optimisation of success probabilities of A. Here 
we will not explicitly use this decomposition further, 
but the constraint will be implemented implicitly by re- 
quiring the n — 1-mode extension to be unitary. 

In our specific 3-mode extension we choose Ai 3 and 
A2.3 such that the new row vectors are orthogonal. By 
multiplying them by the root of the inverse of their re- 
spective norms, \A[ \ and | |, they will be normalised. 
Finding a third orthogonal vector to fill the unitary 
matrix can be done with the complex cross product 
{A[ X A2)*, or in general by choosing a vector at ran- 
dom and orthogonalising it with respect to the given 
ones. 

The dependence of the success probability on the exten- 
sion is 

-2 



Ps = {\A[\\A'2\ 
Therefore, the objective is to 
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minimise/ I a; plA^P (12) 
subject toA'i(A'2)^ = AiAl + Ai^sAl^ = 0. (13) 



The first observation is that the row-scaling by x is al- 
ready included in the norm of the row vectors, leav- 
ing us with one parameter less. By using the phase 
of y, we can assure that AiA^ is real and positive 
and also arg(74i,3) — arg(A2.3) = ±7r. This con- 
strained minimisation problem in Ai 2 and Ai 3 can 
indeed be solved (by using Lagrange's multiplier rule 
and showing constraint qualification) and we find \y\ ~ 



(2(1 - cosip)) 
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Then an optimal solution (phases 



chosen conveniently) is 
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FIG. 2. Basic spatial modes based set-up obtained from translating 
an arbitrary 2x2 core into the language of linear optics. The core 
extension is provided by mixing with a vacuum mode on the cen- 
tral beam splitter. This mode, in turn, has to be post-selected in the 
vacuum state afterwards. The upper and lower beam splitters imple- 
ment the appropriate compensation by damping the by-passed modes 
(which is the same for both modes for the optimal solution which we 
consider). 

The labels at the beam splitters will be used to identify them with the 
respective optical elements in Figs. [3}{5] In general, the parameters 
(i.e., reflectivity and phases) of these elements depend on the gate's 
phase tp. Further, the notation |o)^ and is used for the logical 
0/1-modes of the z-th qubit to avoid confusion with Fock states. 



FIG. 3. Setup for a controlled phase gate on two polarisation encoded 
dual-rail qubits. The logical modes of the two qubits are separated 
and re-united by means of polarising beam splitters (PBS). Replac- 
ing the left and right beam splitters in Fig.|2]by wave plates Ai and 
A2, they become easier to tune to different and provide better sta- 
bility. The lower beam splitter, 7, implements compensation of both, 
|o)-^ and 10)2, modes. \v is taken care of by the partially polarising 
beam splitter (PPBS, allowing for different reflectivities for the two 
polarisations, further explained in Fig.|4ll in the centre. Additionally, 
one of the qubits has to be flipped prior to and after the circuit, which 
here is done by acting with a wave plate on the second qubit. 



with the probabihty of success given by 

ft(^) = (l + |2/P + |y|V2HyF)"' 
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The reflectivities of the compensating beam splitters in 
the by-passed modes have to be chosen such that the 
success probability is constant for all dual-rail states, 
i.e., 



1/4 



(16) 



The success probability Ps{^p) of this gate is shown in 
Fig. [1] Interestingly — and quite surprisingly — the worst suc- 
cess probability is not achieved for the sign-flip {ip = tt), 
but for Lp « 2.05. This means, gates delivering a phase 
shift slightly smaller than tt and thereby generating less en- 
tanglement will not give rise to a larger, but to a smaller suc- 
cess probability. As expected, the success probability for very 
small phases increases and reaches unity for ip = Q: one can 
always do nothing at all with unit success probability. 



C. Integrated quantum photonics realizations 

Sophisticated circuits such as the one shown in Fig.|5]can be 
built from bulk optical elements (mirrors, beamsplitters, etc.) 
Such circuits often require implementation of Sagnac inter- 
ferometers (e.g., Ref. II18I1 ). partia lly p olarizing beam splitters 
(PPBSs) |@], or beam displacers ll9l Isoll to achieve interfer- 
ometric stability. For the most complicated circuits a combi- 
nation of these elements is required. Indeed a (non-optimal) 



implementation of a two-qubit controlled unitary gate used 
a combination of beam displacers and PPBSs However, 
such circuits are extremely challenging to align, are limited 
in performance by the quality of that alignment, and are ulti- 
mately not scalable. 

An alternative approach based on lithographically fabri- 
cated integrated waveguides on a chip has recently been devel- 
oped ilt]. This approach has demonstrated better performance, 
in terms of alignment and stability, as well as miniaturization 
and scalability. The monolithic nature of these devices en- 
ables interferometers to be fabricated with precise phase and 
stability, making the Mach-Zehnder interferometer shown in 
Fig. |2] directly implementable without the need to stabilize 
the optical phase (either actively, or using the Sagnac-type ar- 
chitecture of Fig. |5]l, greatly simplifying the task of making 
complicated circuits: Essentially the circuit one draws on the 
blackboard can be directly 'written' into the circuit. Indeed, 
integrated photonics circuits have been used to implement a 
circuit of several logic gates on four photons, to implement 
a compiled version of Shor's quantum factoring algorithm in 
this way ll20ll . In fact laser direct write techniques have been 
used to 'write' circuits in an even more direct fashion than the 
lithographic approach. 

Another key advantage of integrated quantum photonics for 
the circuits described here is that on-chip phase control can 
be directly integrated with the circuit ll2lll . which could allow 
measurment of the success probability curve (Fig.[T]or Fig.|6] 
for example) to be directly mapped by sweeping the applied 
voltage. 



D. Experimental issues in free space 

In order to render the proposed gates experimentally more 
feasible in free space, some simplifications have to be done 
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work into a single interferometer, is shown in Fig. [3] 

Due to the asymmetry of the core, still one PPBS is used, 
the reflectivity of one polarisation component of which actu- 
ally depends on ip, the other one being 1. Fig. |4] shows how 
a tunable PPBS can be constructed, introducing another inter- 
ferometer. 

Iterating the ideas that led to the compact PPBS implemen- 
tation once more yields a collapsed form of the phase gate 
based on only two PBS and a couple of wave plates. Due to 
the paths of light being very similar, this set-up should also 
be more robust. This is illustrated by Fig. [5] since it has an 
overall Mach-Zehnder interferometer structure to it, by feed- 
ing one displaced sagnac loop into another 



FIG. 4. From left to right: (i) A PPBS implementing a beam splitter 
with polarisation dependent reflectivity, (ii) It is equivalent to an 
interferometer between two PBS where the PPBS' reflectivities are 
incorporated by means of wave plates, Ai and A2. (iii) By identifying 
the two PBS, the interferometer collapses into a closed loop (which 
is more compact and more robust in experimental implementations), 
leaving only one PBS. In the second and third circuit, a polarisation 
flip of the second qubit before and after the circuit is added by using 
a wave plate. 



- tailored to the specific physical implementation at hand. 
Waveguide based setups would not need further simplifica- 
tions because stablity of the interferometers would be ensured 
by the rigid substrate. Implementations more suitable for 
beam displacer based setups are known ITJ. Influenced by the 
gate model, those gates include a controlled 7r-phase gate at 
the core. Operating at lower success probabilities, especially 
at phases (p approaching 0, the probability of success does not 
converge to 1. 

In the following we will discuss simplifications which shall 
allow for easier free-space implementation of the networks in- 
troduced above while still preserving optimality with respect 
to ps- A straightforward set-up on dual -rail encoding that re- 
alises a three-mode unitary and compensates the amplitudes 
in the remaining modes is shown in Fig. [2] It includes one 
interferometer, but the whole gate would sit inside a double 
interferometer, because local unitaries on the input and output 
qubits would require classical interference. Thus, the com- 
plexity of this gate is best described as a nested three-fold in- 
terferometer. In this first stage, the parameters (reflectivities 
and phases) of all five beam splitters depend on ip. 

To get rid of some of the interferometers, polarisation en- 
coding is convenient. Two modes can be united in one spatial 
mode, resulting in inherent stability (neglecting birefringence 
of the optical medium) of some interferometers. Rotations on 
these modes can be carried out easily using wave plates. Be- 
cause the core acts on modes coming from different dual-rail 
qubits, they have to be combined into a single spatial mode 
before. This is achieved with a PBS, thus permuting H and 
V modes. Damping of both by-passed modes can be done si- 
multaneously by a single polarisation insensitive beam split- 
ter coupled to the vacuum. A straightforward translation of 
Fig. |2] into polarisation encoding, thereby collapsing the net- 



E. Simple decomposition 

Here we want to obtain a simple decomposition and explicit 
(^-dependence of the involved elements of the full 3x3 uni- 
tary transformation U e SU{3). To do so, we interpret the 
effective core, acting on the and |i)2 modes in Fig. |2] 
in a different way: in between U\-^ (the unitary beam splitter 
matrix of Ai) and Ux^, there acts a diagonal mode transforma- 
tion such that the first mode is unaffected and the amplitude 
of the second one is damped (due to Ux^, coupling it with re- 
flectivity rv to a vacuum mode which will be projected onto 
the vacuum). 

Now considering the singular value decomposition (SVD) 
of the matrix representing the core transformation, A ^ V ■ 
S • lyt we can identify V = Ux^, W^f = Ux, and S = 
diag{l,ry}. As we have seen earlier, optimal extensions of 
2x2 cores only require global rescaling, which commutes 
with the unitaries involved. Therefore we can use the much 
simpler original form of A in Eqn. and we choose x ~ 1 
and y = — ^/c'-f — 1. 

We find the singular values of A 

cr± = y/l + 2 sin I ± 2(2 - 2 COS ip)^/* cos (17) 

Global rescaling amounts to fixing the largest singular value 
to unity, so the new singular values are 1 and ry = (7-/(7+. 

Due to dct Uxi = det Ux2 1 we need to attach a phase to 
one singular value as well in order to apply the identification 
of the matrices introduced above. Then the SVD of A yields 



and ?7a2 = U^^^ ■ e">'+''' with 



-1 1 
-1 -1 



(18) 



arccot 



cot^±((2-2cos^)i/4sin^ 



(19) 

where the order of rows and columns in the matrices is as in 
Fig.|2]from top to bottom. 

The "complex singular values" are 1 (the first mode is 
not affected) and (T_ /(7-|_ cxp i0+ + (/)_. The latter can be 
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achieved by using the aforementioned coupling to the vac- 
uum with a reflectivity of ry and a phase upon reflection of 
(f)+ + 0- . The further ingredients are the beam spHtters re- 
quired for the "damping" of the by-standers as discussed ear- 
lier By confirming l/a^ = ps in the range < < tt, the 
optimality of this construction is assured. 



III. EVENT-READY GATES 

Coming from post-selected gates, the next step towards 
scalable quantum computation would be to build gates not re- 
quiring measurements on the output modes. Intuitively it is 
clear that the construction of a controlled phase gate in this 
class will be more demanding with respect to the resources 
(such as the number of auxiliary modes and photons, size and 
complexity of the network) involved. 

Especially the number of additional photons will change 
drastically: having had none in the post-selected case of a 
controlled tt -phase gate, two are required in the class of event- 
ready gates. We will use this example as a motivation for a 
detour to discussing a number of different methods that could 
be useful for handling linear optics state preparation. 

To do so, we notice that a controlled 7r-phase gate is more 
constrained than a device that creates EPR pairs from single 
photons. This is meant in the sense that it not only amounts to 
a state transformation from two single photons to an EPR pair, 
but a full unitary transformation on the entire computational 
state space in dual-rail encoding (and creating an EPR pair 
when applied to a certain product input corresponding to the 
product state of two photons). In the following two sections 
we will be concerned with different methods to describe linear 
optics state preparation and will apply them to the specific 
example at hand (i.e., heralded dual-rail EPR pair generations 
from single photons). 

It will turn out, that the construction of an EPR pair out of 
single photons by means of linear optics, vacuum modes, one 
additional photon, and detectors is not possible. Of course, 
directly solving the polynomial equations in the matrix ele- 
ments of A (generalisation of Eqns. (O to Q) will yield the 
same result - no solutions unless two additional photons are 
involved. Having excluded the cases of zero and one auxiliary 
photons, a set-up with two of them is possible, proven by the 
existence of such a scheme (EPR construction ll22ll as well as 
controlled-Z gate iH]). 



A. State transformations 

An obvious way of looking at states of exactly 2 photons in 
TO bosonic field modes is the following. Such a state vector 
can be written as 

IV'm) = P(a^)|vac) 

rn 

= Mijaja]|vac) = (a^)^MaVac), (20) 



where M is a symmetric m xm matrix [? ]. The application of 
a unitary mode-transformation U — representing a linear opti- 
cal network — is reflected by 

\^m) ^ \^M') = (C/at)^7\f(C/at)|vac) (21) 
= (a^)^A/'at|vac) (22) 

with M' = U'^MU clearly again being symmetric. As a 
special case of the singular value decomposition ll23ll . a diag- 
onal M' can be achieved, given an arbitrary input state vector 

Now let us choose U such that M' is diagonal. Then, la- 
belled by 

1^' = rank(M'), (23) 

there are m different classes of states lfl6l [itIi . in each of 
which is the states are composed by superpositions of 2 pho- 
tons in either of v' modes. These classes are separated by lin- 
ear optical mode transformations requiring additional modes. 
Decreasing the rank is possible by allowing for auxiliary vac- 
uum modes. However, to increase the rank by 1 one additional 
photon is required. 

Further, the state matrix M of two single photons on four 
modes has rank = 2 while an EPR pair corresponds to a 
matrix with rank v ^ A. Therefore, the desired state transfor- 
mation requires at least two additional single photons. 

B. Polynomial factorisation 

An alternative approach is the following iflill : The polyno- 
mial describing the objective state vector = P(a^^)|vac) 
with 

P(at) = 2-^1^ {a\al + 44) (24) 

does not factorise over (D. Using Lemmata 1 and 2 from 
Ref. ll24ll . the property of factorisation of a bivariate polyno- 
mial 

m 

P{x,y) = ^ Pijx'-y' (25) 

over (D can be tested by assessing the rank of a complex 
2to(2to — 1) X (to + _i){2m — 1) matrix. Further, apply- 
ing Lemma 7 of Ref. lIZSIl , this technique can be extended to 
multi-variate polynomials. Now, a state can be constructed 
from a product state using linear optical gate arrays iff the cor- 
responding polynomial is factorisable. In the case mentioned 
before (dual-rail EPR pair, so four variables), one can use the 
resulting 12x9 matrix to confirm in the language of polyno- 
mials of creation operators that additional resources are in fact 
required. 

IV. TOFFOLI GATES 

In the same way as above, we can consider a generalised 
Toffoli gate, the effect of which on the computational basis 
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FIG. 5. Compact implementation of a controlled phase gate by using 
a single loop to implement the central PPBS and the compensation 
beam splitters simultaneously. Additionally, the two PBS are iden- 
tified, resulting in a second loop. All omitted modes are initialised 
in the vacuum and post-selected in the vacuum state (which will be 
achieved in practice by counting the photons in the other output). 



0.15 




0.05 



FIG. 6. Optimal success probabilities of generalised Toffoli gates. 
The features exhibited by ps (ip) are similar to the ones observed at 
the controlled phase gate (Fig.[T]l: There is a shallow dip between 
ip = n/2 and ip = n below Psiir) = Pb(ti/2) and a steep in- 
cline (more pronounced than for the controlled phase gate) for small 
phases towards ps(0) = 1. 



realized as dual-rail encodings can be described by the unitary 

[/ = diag(l,l,l,l,l,l,l,e*^). (26) 

The solutions to the polynomial equations describing the ac- 
tion on the three-mode core — up to mode permutations — can 
be parametrised by x, ?/ e ]R and are given by the matrix 





' 1 












xy 
1 


y 




X 





1 



(27) 



This has to be understood similar to the 2-mode core used by 
the controlled phase gate. However, here we do not solve the 
full optimisation problem, but only consider global rescaling. 

For a unitary extension all singular values of A have to be 
at most 1. In order to avoid to formulate cubic singular values 



explicitly, we use the following constraints. Let paa-> W = 
det(AA^ — AI3) be the characteristic polynomial of AA^, the 
roots Ai,2,3 of which are the squared singular values of A. By 
requiring PAAt(l) = 0, one of the singular values has to be 
1. For the roots of p^At are real-valued and non-negative, the 
condition that all other singular values are not larger than 1 is 
equivalent with the condition that all derivatives of p^^t have 
the same sign at A = 1. More formally, this results in further 
constraints of the form 



dA 



> 



(28) 



for 1 < fc < n where a = 0, 1 is fixed by the condition 
of = n. For tp = tt the optimal ps compliant with these 
conditions is 



ft(7r) = 1 + 3 (2I/3 - 22/3) « 1/57. 



(29) 



See Fig. |6] for the maximum success probability in the range 
Q < Lp < TT. The corresponding networks could be con- 
structed in the same way as above. However, they would 
consist of 3-mode cores (an interferometer composed of three 
partially polarising beam splitters) inside separate interferom- 
eters for each of the 3 qubits. For free space experiments more 
appealing approaches for the specific choice of </? = tt are pre- 
sented in Refs. ll26ll27ll . but leading only to success probabil- 
ities of at most 1/72. 



V. REMARKS ON PROCESS TOMOGRAPHY 

Process tomography amounts to characterizing (or some- 
times certifying) an unknown physical process. In practice, 
the task is to identify that completely positive map that is 
closest to the date with respect to some meaningful figure of 
merit. To accomplish this task, one has to consider a tomo- 
graphically complete set of inputs and look at Hilbert-Schmidt 
scalar products of the output with observables, to faithfully re- 
construct the matrix form of the channel Isol [31I1 . Practically, 
the closest physical process can then be found by solving a 
convex optimization problem. Formally equivalently, and in 
instances in an experimentally simpler fashion, one can, in- 
stead of sending in a full set of input states, submit half of 
a single fixed maximally entangled state, and hence recon- 
struct the channel from the Choi matrix, then referrred to as 
entanglement-assisted process tomography. 

The latter technique can clearly also be applied in case of a 
postselected quantum gate like a phase gate. Yet, even without 
entangled inputs and including the actual measurement, one 
can reconstruct the resulting POVM elements, to which es- 
sentially postselected gates amount to when one faithfully in- 
cludes also the actual measurement in the black-box descrip- 
tion of the process. If one has a well-characterized source at 
hand, then the statistics of 



P3,k = tr[pjAfc], 



(30) 



uniquely characterize the process, where {pj} form a to- 
mographically complete well-characterized input set, and 
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Ai,. .., Ak > constitute a POVM, i.e., 

K 

Y.Ak = l (31) 

fc=i 

(for an experimental realization of such an approach, see Ref. 
113211 ). In practice, one uses methods of convex optimization to 
identify the closest physical process to the given data. 

For the purposes of the present work, one of the outcomes 
fc = 1, that is, a specific pattern Ai of detection, then gives 
rise to the actual postselected linear optical quantum gate. 
In this way, one can reconstruct postselected quantum gates, 
without using an ancilla-based approach, even faithfully in- 
cluding the final measurement as part of the process. 



VI. CONCLUSION 

We have shown how to obtain the maximum probability 
of success of controlled phase- and Toffoli-like gates in the 



class of post-selected linear optics dual-rail gates without ad- 
ditional photons. Further, constructions of networks for the 
smaller gates suitable for experimental implementation have 
been given, and techniques elaborated upon that allow for the 
assessment of the possibility of certain linear optical schemes. 
For further progress concerning the eventual full optimization 
of linear optical processes it would be interesting to investi- 
gate (i) optimal constructions with respect to a class of gates 
inherently incorporating such experimental constraints or (ii) 
identify further decomposition techniques from a given linear 
optics mode transformation into suitable physical networks, 
respecting these constraints. 
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